The Electron Spectral Function in Two-Dimensional Fractionalized Phases 
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We study the electron spectral function of various zero-temperature spin-charge separated phases 
in two dimensions. In these phases, the electron is not a fundamental excitation of the system, but 
rather "decays" into a spin-1/2 chargeless fermion (the spinon) and a spinless charge e boson (the 
chargon). Using low-energy effective theories for the spinous (d-wave pairing plus possible Neel 
order), and the chargons (condensed or quantum disordered bosons), we explore three phases of 
possible relevance to the cuprate superconductors: 1) AF* , a fractionalized antiferromagnet where 
the spinous are paired into a state with long-ranged Neel order and the chargons are 1/2-filled 
and (Mott) insulating, 2) the nodal liquid, a fractionalized insulator where the spinous are d-wave 
paired and the chargons are uncondensed, and 3) the d-wave superconductor, where the chargons 
are condensed and the spinous retain a d-wave gap. Working within the Z2 gauge theory of such 
fractionalized phases, our results should be valid at scales below the vison gap. However, on a 
phenomenological level, our results should apply to any spin-charge separated system where the 
excitations have these low-energy effective forms. Comparison with ARPES data in the undoped, 
pseudogapped, and superconducting regions is made. 



I. INTRODUCTION 

Ideas of spin-charge separation have long been consid- 
ered in relation to the cuprate high-Tc materials following 
P.W. Anderson's original suggestions Phenomenolog- 
ically, the assumption that the electron "breaks apart" 
leads to fairly simple explanations for some otherwise 
puzzling aspects of these materials. Attempts to formu- 
late this rather elegant idea into a well-defined theory of 
electrons living in two or more spatial dimensions have 
historically been plagued with problems. A recently in- 
troduced Z2 gauge theory of strongly correlated electron 
systems indeed contains both spin-charge separated 
and spin-charge confined phases, and we work here within 
this formulation. 

Among the host of puzzling experimental properties 
of these materials, we wish to concentrate here on an- 
gle resolved photoemission spectroscopy (ARPES) exper- 
iments which in recent years have reached an unprece- 
dented level of resolution. With this increased clarity of 
data has come increased confusion in theoretical interpre- 
tation. In particular, it seems quite difficult to explain 
the ARPES lineshape in the pseudogap regime within 
Fermi liquid theory. In fact, any conventional quasipar- 
ticle description would seem to predict a sharp peak in 
the spectral function, A{k,bj), at ijj{k) for some k in the 
Brillouin zone. The data in the underdoped compounds 
in their non-superconducting state, on the other hand, 
show only broad and sometimes step-like features. In- 
creased energy and momentum resolution has made the 
contrast with the superconducting state, where a sharp 
peak does emerge, more striking and has led to further 
doubts about the quasiparticle description of the pseudo- 
gap state. As argued elsewhere j^], this contrast between 
the pseudogap and superconducting lineshapes suggests 



that the pseudogap region could be dominated by a zero- 
temperature fractionalized phase. In addition, recent re- 
sults in the superconducting state suggest a connection 
between the weight under the superconducting quasipar- 
ticle peak and the condensate density Q. This result 
seems rather mysterious from a Fermi-liquid point of 
view, but, as we later show, may have a simple expla- 
nation in terms of separated spin and charge degrees of 
freedom. ARPES experiments on undoped compounds 
also show broad spectral features rather than well-defined 
quasiparticle peaks, which has led us to consider the pos- 
sibility of a fractionalized antiferromagnet, dubbed AF* . 
However, the spectral function does show signs of "sharp- 
ening up" as the system is overdoped, suggesting that 
there may be a quantum confinement critical point in 
the cuprate phase diagram, as shown in Figure |l|. 
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FIG. 1. Schematic phase diagram for the high Tc cuprates. 



We wish here to explore in more detail the conse- 
quences of these spin-charge separation ideas for the 
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single-electron spectral function of the cuprate materi- 
als at low doping. Working with a fairly simple theory of 
low-energy spin and charge excitations in a fractionalized 
phase, we will find qualitative agreement with ARPES 
data in the pseudogap and superconducting phases, as 
well as in the undoped insulator. Although the theory 
used here has been analyzed and motivated from a variety 
of standpoints elsewhere ||^,^ , we hope to make clear its 
reasonableness on purely phcnomcnological grounds. We 
begin, then, from a zero-temperature theory of d-wave 
paired spinous and charge e bosons. The bosons can 
have a zero-temperature phase transition between con- 
densed and quantum disordered phases. We explore here 
quantitatively the single-electron spectral function in the 
a; = spin-charge separated antiferromagnet {AF*), the 
nodal liquid (to be identified with the pseudogap phase) , 
and the superconductor. 



II. THE MODEL 

We briefly recapitulate the phase diagram of the 
cuprates in terms of the Z2 gauge theory introduced else- 
where The theory contains spinon and chargon de- 
grees of freedom, coupled to a Z2 gauge field in two spa- 
tial dimensions. We begin with the square lattice Hamil- 
tonian: 

i i 

-^E-s-^En^^' (1) 

<ij> □ □ 

where the electron operator is a product of spinon and 
chargon operators: Cia = bifia- The term with coupling 
K is allowed by symmetry and can arise from integrating 
out the very high energy chargons, making this an effec- 
tive theory of the low energy charge degrees of freedom. 
The spinon pairing Ay is taken to be d-wave: 

^ \ -A along y, ^ ' 

and the spin operator is = J2k fl+n'^fk- ^ is the 
mean-field Neel order parameter and is non-zero only 
within the antiferromagnetic phase. The U term is a 
Hubbard-like interaction for (1 — x) chargons per unit 
cell. At zero temperature and as a function of K/h, the 
gauge field has a transition between confining and de- 
confining phases Deep within the deconfining phase, 
we may set a^j = 1 on all links and we are left with 
decoupled spinons and chargons: 



<v> 

+Uj2[n,-il-x)r+Y,NJr- (3) 

i i 

Fluctuations of ct^ can be taken into account by con- 
sidering vortices in the Ising gauge field which have been 
dubbed "visons" . (A plaquette which contains a vison 
has On ^fj = ^1-) The deconfining phase of the Z2 
gauge field is characterized by a gap to these vison exci- 
tations and, as we see above, the electron degrees of free- 
dom are fractionalized in this phase. The zero tempera- 
ture confining phase of the Z2 gauge field is a condensate 
of these vison excitations and "glues together" spinons 
and chargons to form electrons. A "quantum confine- 
ment critical point" separates these two zero-temperature 
phases, as discussed elsewhere [1). At finite temperatures 
above the fractionalized zero-temperature phase, we ex- 
pect vison excitations to exist in two dimensions, leading 
to interactions between the chargons and spinons. How- 
ever, at temperatures much smaller than the vison gap 
rpvison (^^ff\y{Q\y jjj simplest theories of the quantum 
critical point would be of the same order as the pseu- 
dogap temperature, T* ), we expect the low energy de- 
grees of freedom to be those of the fractionalized phase: 
spinons and chargons weakly interacting through Doppler 
shift terms, which we ignore. 

We briefly discuss the phases shown pictorially in Fig- 
ure |l|. In the Xc < X < xqcp range, starting at tempera- 
tures much less than the vison gap and lowering the tem- 
perature, the bosonic chargons should go from a phase 
where they are phase incoherent to one where they are 
condensed. Below the chargon condensation tempera- 
ture, the system is superconducting; this is Tc. Through- 
out, the spinons maintain a d-wave pairing (presumably 
due to antiferromagnetic interactions of strength J) and 
experience no phase transition, but rather a crossover 
at their pairing scale, T*. Starting instead at zero- 
temperature and zero doping, we are in a spin-charge 
separated phase which is also an antiferromagnetic Mott 
insulator with long-range Neel order. Upon increasing 
the doping, staying at zero-temperature, we presumably 
enter a complicated charge-ordered insulating state of the 
chargons and destroy the long-range Neel order of the 
spinons. This is the zero-temperature phase believed to 
dominate the pseudogap region. We expect impurities 
and thermal fluctuations to destroy static charge order, 
but inhomogeneous effects could still be an important 
high-energy presence, leading to e.g, stripes. As the dop- 
ing is further increased, the chargons presumably con- 
dense at zero-temperature into a superconducting state. 
After the destruction of Neel order, the spinons are qual- 
itatively the same in this doping range and maintain a d- 
wave gap of order T* . Throughout this zero-temperature 
region, the chargons and spinons are decoupled (since we 
are to the left of xqcp)- At a; = xqcp, the Ising gauge 
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field becomes confining and the chargons and spinons are 
bound together to form electrons, presumably in a Fermi 
liquid phase. 

We turn our attention now to the spectral function 
defined in terms of the electron Green function: 



A{k,Lu) = -—ImG{k,uj) 



(4) 



Since at temperatures well below the vison gap, we ex- 
pect a description of the system in terms of free chargons 
and spinons to capture the low-energy physics, we use the 
Hamiltonian in Equation |^, which is a sum of spinon and 
chargon Hamiltonians: H{c\c) ~ Hb{b\b) + Hf{f^ ,f). 
Within this construction, it is possible to write the elec- 
tron Green function as a product of chargon and spinon 
Green functions: 



G(r, 



= (T,c(r,r)c^(O,0)), (5) 
= (T.6(r,r)&t(0,0))(r./(r,r)/t(0,0)}, (6) 
= Gb{r,T)Gfir,T), (7) 



with T = it, the imaginary time. The problem of cal- 
culating the spectral lincshape in spin-charge separated 
phases now becomes one of calculating the spinon and 
chargon Green functions. We consider these two degrees 
of freedom in turn, discussing values of various parame- 
ters in each phase. 



A. Spinons 

We briefly describe the phases of the spinon model. 
Consider first the N = phase which describes spinons 
with a d-wave paring amplitude, Afc. In this spin-charge 
separated construction, superconductivity is dependent 
only on the charge degrees of freedom. When the bosonic 
chargons are condensed ((6) 7^ 0), we are in a BCS d-wave 
superconductor and the spinons are simply neutralized 
BCS quasiparticles. When the chargons lack phase co- 
herence, we are in a phase with no superconductivity, but 
with a d-wave gap to any excitation with spin 1/2, called 
elsewhere the nodal liquid When ^ ^ 0, spinon- 
antispinon pairs condense forming a state with long-range 
antiferromagnetic order, but still containing free spinon 
excitations above a gap (of order J) which are separated 
from the chargons due to the vison gap. This spin-charge 
separated antiferromagnet has been dubbed AF* |^ . 

The spinon piece of the Hamiltonian in Equation ^ is 
quadratic in the spinon operators, and we may diagonal- 
ize it using a Bogoliubov-type transformation. Setting 
gN = N^z and working in units of the lattice constant, 
we obtain: 



(8) 
(9) 



with. 



Efc = —ts{coskx + cosky), 
Afc — A{coskx — cosky), 



ak,a = Ukdk^a + OiVkd^_k_^, 

i cos6'fc ; v1^\ - i cos6'fc. 



COSffc = 



where 



/l2 - 1 



dk,a — Akfk.a + 0:Bkfk+n,a, 



COS (pk 



BI = ^ — COS( 



COS (pk = 



Ek 



(10) 
(11) 



(12) 
(13) 
(14) 



(15) 
(16) 

(17) 



is a Hartree-Fock-type spin density wave operator at mo- 
mentum TT appropriate to commensurate antiferromag- 
netic order with sublattice magnetization A^o Note 
that when A^o = 0, the Hamiltonian for is the same as 
the effective BCS Hamiltonian for a d-wave superconduc- 
tor. Indeed, when the chargons condense, these become 
the Bogoliubov d-wave quasiparticles. 

At zero temperature, we have for the spinon correla- 
tion function: 



0a,l3 

2 



1 



Ek 



(18) 



We see that the spinon spectrum now has a gap of A^o at 



kr 



^, as we would expect in the Neel state, as 



well as a d-wave gap whose maximum is 2A . We expect 
this spinon theory to work qualitatively at all temper- 
atures well below T* < T^^''°^, where the spinons are 
strongly paired and the low energy degrees of freedom 
are fractionalized. In particular, one should note that in 
the absence of chargon-spinon interactions, the spinons 
do not notice Tc- 

The parameters t^ and A can be set by the experimen- 
tally determined ratio: 



k 

A 



VA ' 



(19) 



which ranges from - 14 in YBCO to ~ 20 in BSCCO 
near optimal dopings 0. We expect this ratio to be of 
this order throughout the pseudogap phase. At zero dop- 
ing, I TV I is on the order of J. 



B. Chargons 

Once liberated from their fermionic statistics, the 
charge degrees of freedom behave as bosons of charge 
e hopping on a 2d square lattice, as in Equation ^. At 
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half-filling and zero temperature, we expect that in the 
limit U/tc » 1, the system forms a Mott insulator, 
while in the limit tc/U » 1, the bosons form a su- 
perfluid. This can be described by the 2-|-l-dimensional 
quantum XY model which has two phases, a supercon- 
ducting phase and a quantum disordered, Mott insulating 
phase. Being concerned primarily with "normal state" 
(i.e. non-superconducting) properties, consider the in- 
sulating phase where U/tc » 1- Excitations of this 
phase are doubly-occupied sites which are "massive" (i.e. 
gapped) and may propagate. These excitations as well 
as the excitations within the superfluid phase are well- 
described by the soft-spin continuum Landau-Ginzburg 
action, replacing the chargon operator with a complex 
field. 

k ^ b{r), (20) 

= i \drb\^ + 4|V6|2 + H|5|2 + ^(|5|2)2^ (21) 

When > 0, the bosons are quantum disordered and the 
chargon system is insulating. When /i < , the chargons 
condense, forming a superconductor with \ {b)\'^ = fi/Au — 
no, where uq is the condensate density. Fluctuations 
around this new minimum are described (to quadratic 
order) by the action: 

6(r,t) = (6)+6(r,i), (22) 

6 = &l+^&2, (23) 
A = i (^A)' + ^(V&i)' + ^(^i)' 

+ ^ (9,62)2 + ^(V62)2, (24) 

with, = -2fi. (25) 

Starting in the superconducting phase and increasing the 
chemical potential toward /i = 0, the order parameter 
(and therefore the condensate density) vanishes at the 
transition. 

Away from half-filling in the presence of long-range 
Coulomb interactions or disorder, we expect the unoc- 
cupied sites to form some crystal. Even in the case of 
zero disorder, the underlying lattice makes characteriza- 
tion of this phase difficult. One way to gain intuition for 
this regime is by reformulating the problem in terms of 
vortices in the chargon phase On physical grounds, 
we expect that the strong coupling of the charge degrees 
of freedom will lead to complicated charge ordered states 
at zero temperature when the number of bosons is in- 
commensurate with the underlying lattice, and that with 
increasing doping, the system should eventually pass into 
a zero-temperature superconducting state. The location 
of the transition, Xc, and the nature of the exact ground 
state for x < Xc will depend sensitively on the chargon 
interactions and lattice commensurability effects. Lack- 
ing a more detailed theory of chargon solidification away 
from half filling in a fractionalized phase, we will use 
the XY model defined above to describe the low-energy 



degrees of freedom at low temperatures in the fraction- 
alized phases. Our main motivation is simplicity: the 
2+ld XY model contains both a quantum disordered and 
a superconducting phase of bosons, as the more correct 
theory of the x > boson system should. Although this 
description is obviously inadequate to describe the zero- 
temperature phases away from half-filling as well as the 
detailed critical properties of the transition, we note that 
for a perfectly clean system, the physics at length scales 
shorter than the mean hole spacing should be those of 
the half-filled system. At dopings of, say, 5 percent, this 
length is about 5 lattice spacings. In the correspond- 
ing energy range, the 2-1-1 XY model should capture the 
correct physics. We note that ARPES is an intermedi- 
ate energy probe, although the energies corresponding to 
moderate dopings may still be too high. Since we are con- 
cerned here with general features of spectral function in 
each phase, we work with this phenomenological Landau- 
Ginzburg model, hoping to capture the correct physics. 

Therefore, in the charge disordered (nodal liquid) 
phase at temperatures much smaller than the vison gap, 
we use Equation |2^ and find for the chargon correlation 
function (setting ^ = w? > and ignoring u to lowest 
order): 

(bibk) = — , (26) 
Lul^m^ + v^\k\'^. (27) 

We briefly discuss the parameters in the model, m and 
V. In the underdoped regime near the critical point, Xc, 
we expect the chargon gap to be quite small, while at 
half-filling (in the parent insulator) the chargon gap is 
rather large, the charge gap for these materials being 
on the order of electron volts. Very little can be said 
about the velocity v without a more detailed microscopic 
theory. Working in units where fc is a dimensionless 
wavenumber, v is an energy scale, and we take it to be 
moderately larger than the spinon kinetic scale, is, (but 
of the same order) effectively giving the charge excita- 
tions a larger bandwidth than the spinous. 

In the ordered phase (^ < ), the bosons are super- 
conducting and at T = the chargon-chargon correlation 
function has the property: 

, ^'"^ (bHr,t)bir\t'))^\{b)\'=no, (28) 

r — r — > OO 

the condensate density. The quantity (&) is precisely 
the order parameter of the superconductor. At dopings 
x > as temperature increases within the superfluid 
phase, we expect phase fluctuations to reduce this quan- 
tity, eventually causing it to vanish at T = Tc{x). At 
zero temperature within the superfluid phase, as doping 
is decreased it will vanish a,t x = Xc- As discussed above, 
the details of the T = transition will be governed by the 
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universality class of the true doping-dependent chargon 
theory. 

Within the superconducting phase, we may model the 
bosons with Equation |2^, which results in the following 
general form for the chargon spectral function: 

(6t(r,i)6(r',t')) - m{T) + {h\r ,t)h{r\t')) . (29) 

The fluctuations of the chargon field will be dominated 
by the detailed interactions between the chargons. In 
contrast to the Cooper pairs in a standard BCS super- 
conductor, the bosonic chargons should be strongly inter- 
acting, given that their uncondensed phase is controlled 
by Mott-insulating physics. However, we expect that at 
energies larger than the condensation temperature, Tc, 
the chargon fluctuations should be the same as in the 
pseudogap state. 



III. SPECTRAL FUNCTION 

Given the spinon and chargon correlation functions, we 
can compute the electron spectral function, assuming no 
interactions between chargons and spinous, using the re- 
lations in Eqns. ^ and ^. The result at zero temperature 
is: 



A{k,uj) 



+ (4 A) {hi_^bk-^)5{Lo + LOk-^ + E^)\, (30) 
= A+(fc,w) + A_(fc,w). (31) 

Because it measures electrons ejected from the sample, 
the ARPES intensity (up to matrix element effects) mea- 
sures the occupied part of the spectral function. A- {k, lu) 
. At temperatures far below the vison gap, the assump- 
tion of no chargon-spinon interactions should be valid. 
At energies, lo, larger than the temperature, the use of 
zero-temperature results should be valid. Most of the 
ARPES data of interest are done at temperatures below 
100 K, which translates to an energy of 10 meV or less, 
close to the resolution of the instruments and certainly 
smaller than any features in the "normal state" spectra. 
This justifies use of the zero-temperature spectral func- 
tion in our low-energy model. We therefore compare this 
A-{k,uj) with the ARPES data in each of the following 
phases: AF* , nodal liquid (pseudogap), and d-wave su- 
perconductor. Although Equation ^ is quite simple, it 
nevertheless is not analytically integrable for arbitrary k 
and LU. In the following sections, we present the results of 
numerical integrations of this function and plot the resul- 
tant A^{k,uj) at fixed k (EDC) and at fixed u (MDC), 
along the momentum cuts shown in Figure |. For the nu- 
merical integration, we approximate the delta function in 



Eqn. ^ by a Lorentzian of small width (0.0125 eV) for 
the energy distribution curves. This leads to small "tails" 
in these curves at small binding energies (near turn-on). 
Since we would like to explore the momentum distribu- 
tion curves at these small turn-on energies, we need to 
avoid measuring mostly these Lorentzian tails. To this 
end, for the MDCs, we use instead a box-like delta func- 
tion: 



Six) 



1/e for - e/2 < X < e, 
else. 



(32) 



with e = 10 meV. Values of various parameters (such as 
A and v ) will be given in each section. Some "Fermi 
surface" properties of the spinous should be discussed. 
While the pairing terms technically destroy any true 
Fermi surface, the number of spinous at a given momen- 
tum still drops off for cj^ > and is sensitive to the 
minimum of the spinon dispersion, , which we can call 
kf. Along the (tTjTt) direction, this minimum occurs at 
k — (7r/2,7r/2), while along the (tt, 0) direction, the lo- 
cation of the minimum depends on the relative values of 
ts and A and will be discussed in each phase. 




FIG. 2. Momentum cuts used for plots of A-{k, cj), show- 
ing the approximate location of the "Fermi surface" for our 
model. Cut #1 (used for MDCs and EDCs) is along the line 
kx ~ ky near the nodal point, cut #2 (used for EDCs) is along 
ky = near the antinodal point. 



A. (Fractionalized) Antiferromagnet 

In this phase, the spinous are particle-hole paired into 
an antiferromagnet (with single spinous above the gap) 
and the chargons are gapped into a Mott insulating 
phase. However, because these two particles propagate 
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as separate excitations, we expect an electron injected 
into the system to "decay" or fractionalize into these two 
constituents. We therefore expect the spectral function 
at temperatures much lower than the vison gap to be 
broad, without the delta- function peak at some k and w 
which one finds when the underlying phase has electron- 
like elementary excitations. 



(b) Cut #2 




kx=1.8 

kx=2.0 

kx=2.2(=kF) 

kx=2.4 

kx=2.6 



-2 -1.9 -1. 



-1.7 -1.6 -1.5 -1.4 -1.3 -1.2 -1.1 
Binding Energy (eV) 



(c) Cut #1 



For the spinous, we expect that both A and Nq are of 
the order of J ~ ts/2; we take tg ~ O.beV. The chargon 
gap, m, is expected to be fairly large, on the order of 
an electron volt. With this in mind, we plot the elec- 
tron spectral function in Figure || in the AF* phase with 
iVo = A = 0.25 eV, U = 0.5 eV, v = 2.5 eV, and to = 1 
eV. The shapes of the curves are not sensitively depen- 
dent on any of these parameters. For this ratio of tg to 
A, the minimum of the spinon energy, E^, along cvit ^2 
occurs at ~ 2.2. We plot the MDC along cut #1, at 
the energy uj — —1.30 eV. This is slightly larger than the 
minimum binding energy of + to = 1.25 cV. 



(a) Cut #1 



kx=ky=1.4 

kx=ky=1.5 

kx=ky=pi/2 

kx=ky=1.65 

kx=ky=1.75 




-1.9 -1.8 -1.7 -1.6 -1.5 -1.4 -1.3 -1.2 -1.1 
Binding Energy (eV) 




1.2 



1.4 



1.6 
kx=l<y 



1.8 



2.2 



FIG. 3. A-{k,ui) at zero temperature in the AF* phase. 
Plotted are; EDCs along (a) cut #1, and (b) cut #2 and an 
MDC (c) along cut #1 at energy lo = —1.30 eV. The momen- 
tum space cuts are shown in Figure ^. 



A few features of these curves should be pointed out. 
First, all EDCs are quite smeared, with no peaks. This 
indeed mimics one of the features of ARPES data on 
the undoped compounds ]lO|] . Detailed comparison with 
the binding energies in the ARPES data for the undoped 
compound is complicated by the fact that determining 
the "Fermi level" of these compounds is not as straight- 
forward as in the doped materials. In the case of our 
EDC plots, the binding energy should be used to note 
that the leading edge along cut #1 indeed has a smaller 
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gap than that along cut #2. In fact, the difference be- 
tween the two is just the factor of 2 A. That this feature 
of the leading edge will track Aj^, the d-wave gap, can be 
seen from Equation |3^. This is consistent with the ex- 
perimentally determined "remnant Fermi surface" with a 
d-wave character found in the undoped compounds . 
In contrast to the EDCs, the MDC shows a sharp feature. 
The detailed shape of the MDC curves may be influenced 
by the specifics of the model used. 



A_[fc = (7r,0),w ~ A] ~ 

f d\ 1 . ^F[ql-ql] \ ~ 

J (27r)2 Avq y 2Eq j 1' y J 

1 r°° 

-7^ dq5[uJ + vq + 2l^ + 0{q^)] (37) 
^-l-,9i-.-2A). (38) 



B. Pseudogap (Nodal Liquid) 

Because this zero temperature phase is also fraction- 
alized, we expect broad spectral functions in this region. 
The spinous are paired into d-wave singlets, leading to 
the spin-gap. To be precise, we calculate the spectral 
function at zero-temperature using the XY model de- 
scribed earlier. We expect a low-energy theory of quan- 
tum disordered chargons to work qualitatively for the 
entire pseudogap region, provided T « At fi- 

nite temperatures, the zero-temperature spectral func- 
tion can only be expected to capture features at energies 
larger than T. We would therefore like to compare our 
spectral function in this phase with ARPES data in the 
underdoped compounds at T* » T > Tc- The char- 
gons in this region should be dominated by their zero- 
temperature critical point. Here, we use the critical 2-|-ld 
XY theory for the chargons described previously, again 
noting that this will not describe in detail the true finite- 
doping critical point but will hopefully give an adequate 
effective theory for the low-energy excitations. 

As an illustrative calculation, we may analytically per- 
form the convolution integral in Equation |3^ for k = 
(7r/2,7r/2) and k = {tt,0) at small w, exactly at the XY 
critical point, m = 0, of Eqns.E6[]27^. Equation reads: 



(Pq 



1 - 



Cfc- 



k — q 



5{u + Ek-q+UJq). (33) 



At the node, the spinon spectrum may be linearized for 
small momentum and we find (after rotating to momenta 
parallel and perpendicular to the nodal direction and set- 
ting tg = /S. = v for simplicity) : 



A^[k = (7r/2,7r/2),w small] ~ 

d^q 1 A \ xi , - , ^ 



1 

Bttu 



dq5[L0 + [v + v)q\ = 



1 



87rv(w + v) 



(34) 

i-uj). (35) 



At the antinode, the spinon spectral function is quadratic 
above the gap, E^—q 



Eg = 2AJ1 



We see that at these two particular points in fc-space, the 
electron spectral function turns on like a step function, 
not a peak, in qualitative agreement with the ARPES 
results. 

To obtain the electron spectral function at other k and 
uj, we resort to numerical integration. In the underdoped 
region, values of ts/A vary from compound to compound, 
but are of order 10. With this in mind, we set tg = 0.5 
eV, A = 25 meV, v — 2.5 eV. Also, for the purposes 
of performing the integration with the Lorentzian delta- 
function approximation (see discussion at the beginning 
of section III), we regularize the chargon spectrum by 
adding a small mass, m = 12.5 meV, for the energy dis- 
tribution curves shown in Figure |(a) and (b). For the 
momentum distribution curve shown in Figure ^(c) , the 
chargon mass is equal to zero. Again, we find no sensitive 
dependence on the exact values of these parameters. For 
this value of tg/ A, the minimum value of the spinon en- 
ergy, Ej^, along cut #2 occurs at ~ 3.0. For the MDC 
along cut #1, we use a binding energy large enough that 
the width of the approximate delta-function does not in- 
fluence the width of the curve, uj — —40 meV. 



(a) Cut #1 



kx=ky=1 .45 
kx=ky=1 .5 
kx=ky=1 .55 
kx=ky=pi/2 
kx=ky=1 .6 




^g^, and we find: 



-0.3 -0.2 -0.1 0.0 
Binding Energy (eV) 
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(b) Cut #2 



kx=2.3 
kx=2.5 
kx=2.8 
kx=3.0(=kF) 
kx=3.2 




-0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 
Binding Energy (eV) 



(c) Cut #1 



1.2 



1.4 



1.6 
kx=l<y 



1.8 



2.2 



FIG. 4. A-{k,Lj) in the pseudogap phase. Shown are: 
EDCs along (a) cut #1, and (b) cut #2 and an MDC (c) 
along cut #1 at an energy of w = —40 meV. The momentum 
space cuts are shown in Figure H. 



We wish to note the foUowing features of the graphs in 
Figure^. Foremost, the EDCs are indeed quite smeared, 
even near the "Fermi surface" crossings of the spinous, 
but more so in the (tt, 0) direction than in the nodal 
direction, where something peak-ish (though stiU quite 
broad) emerges near (7r/2,7r/2). Also, as we have seen 
analytically, the leading edge in the (tt, 0) direction never 
gets to zero binding energy but instead shows a gap of 
2A = 50meV. It should also be pointed out that as 
one moves along either cut, both sets of EDCs show the 



leading edge moving toward its minimum binding energy 
and then losing weight and/or receding above fc/. Of 
particular interest is the contrast between the EDCs and 
MDC along cut #1 (the nodal direction), where the MDC 
shows a very sharp peak at the node while the EDCs are 
broad and often step-like. The noise at the top of the 
MDC is a consequence of using a "box-like" ^-function 
for this integration. 



C. D-Wave Superconductor 

At low dopings (where T"*^°" >> Tc) when we cool 
below Tc, the bosonic chargons develop phase coherence 
and {b) is non-zero. The single electron correlation func- 
tion in this region then has two pieces, in accordance with 
Eqns. |and|M[ 

G(r,r) = |(6)|2(/(r,r)/t(0,0)) 

+ (6(r,T)5t(O,0))(/(r,T)/t(O,0)), (39) 

giving an occupied portion of the spectral function: 

A^{k,Lo) = no{T){flfk)S{Lu + Ek) 

{flh){bi_^bk-^)5{uj+uu-q + E^). (40) 

Technically, this form is only valid at zero-temperature. 
However, we can see from the electron Green function in 
Eqn. ^ (which is valid at all temperatures much less 
than T"«*°") that throughout the superconducting phase 
at low dopings, we expect a spectral function made up of 
a peak and a background. 

The peak is a product of the condensate density, no (T) , 
and the spinon spectral function. For a given value of 
fc, the peak is located at the BCS quasiparticle energy, 
Ek. Indeed, for non-interacting bosons at zero tempera- 
ture, {b^h) = 0, and we reproduce the BCS quasiparticle 
peak. In contrast with the bosons of BCS theory, we ex- 
pect the chargons to be strongly interacting, leading to a 
non-zero background even at zero temperature | pl| . We 
note that the width of this peak in our simple theory is 
entirely determined by the width of the spinon spectral 
function. Throughout the superconducting state, we ex- 
pect the spinous to act like a two-dimensional Fermi liq- 
uid, leading to a weak temperature-dependent width. 
To the extent that the peak and the background are dis- 
tinguishable objects, the weight under this quasi-particle 
peak should be proportional to the condensate density. 



A^{k,io)=no{T) S{Lu + Ek)^no(T), (41) 



peak 



peak 



and should vanish into the background as T Tc from 
below, without appreciable broadening. 

A comment should be made here regarding the differ- 
ence between the condensate density, no, and the super- 
flxiid stiffness, ps- While for non- interacting bosons these 
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quantities are the same, for interacting bosons they are 
different, even at zero temperature. Besides the effect of 
chargon-chargon interactions on these quantities, there 
is the important effect of the Doppler shift couphng be- 
tween the superfluid and the quasiparticles in the super- 
conducting state. For a d-wave superconductor, the cou- 
phng between quasiparticles and condensate leads to the 
well-known T-linear depletion of the superfluid stiffness 
for small T. The penetration depth, because it measures 
the superfluid stiffness, manifests this dependence near 
T = 0. The condensate density, on the other hand, is 
not directly coupled to the quasiparticles, and therefore 
need not approach T = in the same manner as the 
superfluid stiffness. 



The background in the spectral function comes from 
the second term in Equation |4^, and will be complicated 
by the exact nature of chargon interactions. At ener- 
gies large compared to the condensation temperature (~ 
10 meV), we expect the spectral function to be that of 
the "normal state" above Tc. At low energies, we have 
seen above that there will be a sharp (resolution lim- 
ited) peak in the spectral function, with weight equal to 
the condensate density, located at the spinon gap. It is 
only at intermediate energies, (1-lOmeV), that the de- 
tailed physics of the chargons at their charge-ordering 
critical point becomes important. In the superconductor 
for Tc << T"*''°", we therefore expect a sharp peak whose 
weight is given by the condensate density, superimposed 
on a background which does not change qualitatively as 
one moves from the superconductor to the pseudogap 
phase above Tc. An illustration of the spectral function 
in the superconducting phase is given in Figure |^. 



-0.5 -0.4 -0.3 -0.2 -0.1 0.0 
Binding Energy (eV) 




FIG. 5. For illustrative purposes only, this figure shows a 
Lorentzian peak centered at Ek superimposed on the nodal 
liquid A-{k,uj) for the "Fermi surface" crossing along cut ^2 
from the previous section, k = (3.0, 0). 



IV. CONCLUSIONS 

We have shown here that the following aspects of the 
ARPES data in the cuprate materials can be understood 
by assuming spin-charge separation: 1) the d-wave "pseu- 
dogap" seen above Tc, 2) the lack of sharp quasiparticle 
peaks in the pseudogap phase, 3) the emergence of a very 
sharp quasiparticle peak below Tc, 4) the qualitative tem- 
perature and doping dependence of the weight under this 
quasiparticle peak, as well as the existence within the su- 
perconducting state of a background similar in shape to 
the pseudogap spectra, and 5) the lack of sharp features 
in the undoped parent insulators as well as the d-wave 
character of their "remnant Fermi surface" . We empha- 
size that these results of ARPES in the undoped and 
underdoped compounds are rather hard to account for 
within a conventional picture of quasiparticles with the 
quantum numbers of electrons. 
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